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SUMMARY 

Equations were developed for the torsional motion of a gas-turbine 
engine geared to a helicopter rotor in which the rotor blades were 
hinged to the rotor shaft. The rotor system was simplified to yield 
simple third-order equations that can be used in the analysis of engine 
control. Comparison of the system response calculated from these equa- 
tions with the experimentally observed frequency response of a rotor 
from a 2500-pound helicopter showed satisfactory agreement. 

Calculations showed that the torsional motion arising from the 
hinged construction of the blades contributed to the instability of the 
engine-speed and engine- torque controls. Trends in stability and re- 
sponse of controls with increasing weight of helicopters were 
investigated. 


INTRODUCTION 

Investigation of the control of the speed of a gas turbine geared 
to a propeller is shown in reference 1 to involve the problem of 
controlling the speed of a single rotating body with damping forces 
from the engine and the propeller. When the propeller is replaced by 
a helicopter rotor, however, torsional flexibility between the gas tur- 
bine and the helicopter rotor may be sufficient to establish torsional 
oscillations. In this case, the torsional oscillation may contribute 
an unstabilizing influence to the engine control. 

The torsional flexibility in the helicopter rotor arises from two 
sources: First, the long and highly stressed drive shaft permits con- 

siderable twisting; and second, in many rotors the blades are fastened 
to the hub by lag hinges, which, as shown in figure 1, permit the blades 
to swing relative to the hub in the plane of rotation. Centrifugal 
force tends to hold the blades in a radial position, while aerodynamic 
drag forces tend to make the blades fall behind a radial position. 


2 


NACA TN 3027 


The torsional oscillations of a rotor from a small helicopter (2500 
lb gross wt) were measured and are reported in reference 2. The data 
revealed a badly underdamped resonance at 7.9 radians per second, which 
is a low enough frequency to influence the action of the engine control. 

The investigation reported herein was conducted at the NACA Lewis 
laboratory to determine whether the torsional oscillations in helicopter 
rotors have a detrimental effect on the control of gas-turbine engines 
to which they are geared. The investigation was made by first estab- 
lishing differential equations of motion of the rotor and then verifying 
the equations by comparing the calculated results with the observed os- 
cillations reported in reference 2. 

The effects of helicopter size and rotor weight on rotor dynamics 
were investigated, and the control characteristics of a gas-turbine 
engine in a large helicopter were studied. The control systems con- 
sidered were speed and torque control by fuel flow. 


SYMBOLS 

The following symbols are used in this report: 

radius from center of rotor shaft to drag hinge, ft 

distance from drag hinge to center of gravity of rotor blade, 
ft 

coefficient of characteristic equation 
dimension in geometry of force vectors, ft 
center of gravity 

coefficient in characteristic equation 
dimension in geometry of force vectors, ft 
coefficient of characteristic equation 
dimension in geometry of force vectors, ft 
drag force, lb 
centrifugal force, lb 
tensile force, lb 
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f radius at which F^ acts, ft 

G coefficient of characteristic equation 

hp horsepower of engine 

I e inertia of engine rotor and gearing at rotor shaft speed, 

lb -ft sec 2 

I r polar moment of inertia of rotor blades about center of rotor 

shaft, lb-ft sec 2 

Kq gain of control amplifier, (lb/sec)/(radians/sec) 

Kp loop gain of control 

k a effective spring constant of torque restoring rotor blades to 

radial positions, ft-lb/radian 

kb torsional spring rate of rotor shaft, ft-lb/radian 

kp viscous-drag coefficient of dampers, (ft-lb) (sec)/radian 

ke rate of change of engine torque with engine speed at constant 

fuel flow, (ft-lb) (sec)/radian 

k r rate of change of aerodynamic torque on rotor with change in 

rotor blade speed, (ft-lb) ( sec) /radian 

k s effective spring rate between engine and rotor blades, 

ft-lb/ radian 

ky rate of change of engine torque with fuel flow, ft-lb/(lb 

fuel/sec ) 

kp rate of change of aerodynamic torque with c han ge in rotor 

pitch at constant rotor blade speed, ft-lb/radian 

L lift or gross weight of helicopter, lb 

M mass of rotor blades, (lb)(sec 2 )/ft 

We amplitude of engine speed variation, radians/sec 

W e engine speed, radians/sec 

N r amplitude of rotor speed variation measured at centers of 

gravity of rotor blades, radians/sec 
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N r 


rotor angular velocity measured at centers of gravity of rotor 
blades, radians/sec 



slope of steady-state, speed - fuel-flow characteristic 


P 

% 

Qe 

Q r 

Qs 

r t 

w 

a 

(3 

i 

£sn 


symbol for d/dt 

amplitude of torque force developed by damper, ft-lb 
amplitude of engine-torque variation, ft-lb 
amplitude of rotor-torque variation, ft-lb 
amplitude of shaft-torque variation, ft-lb 
radius to tip of rotor, ft 
amplitude of fuel-flow variation, lb/sec 
amplitude of angular twist of rotor shaft, radians 
amplitude of rotor -pitch variation, radians 

amplitude of angular motion of rotor blades relative to center 
of rotor shaft, radians 

change in steady-state value of drag angle at constant rotor 
speed following change in rotor pitch, radians 


e 


^d 




T r 


amplitude of angular deflection between engine and helicopter 
rotor blades, radians 

time constant of fuel system, sec 

time constant of derivative action of control, sec 
engine time constant, sec 

time constant of integrator action of control, sec 
rotor-blade time constant, I r /-k r , sec 

frequency of sinusoidal disturbance in rotor pitch, radians/sec 
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undamped natural frequency of engine-rotor system, radians/sec 
torsional pendulum frequency of rotor, , radians/sec 



EQUATIONS OF MOTION FOR ROTOR 


In this section of the report, a simplified model of the rotor- 
engine system is developed and the linear differential equations of 
motion are established. The important transfer functions are presented 
and examined. 


A model of the rotor was used to obtain equations simple enough to 
be easily represented by an electronic analog. More complete statements 
of rotor -blade motion during forward flight, including the flexibility 
of the rotor-bearing mounts, appear in standard texts, such as refer- 
ence 3. They are too complex for use in control study but could, per- 
haps, be simplified. 

The articulated rotor system of figure 1 was considered to consist 
two coaxial bodies in pure rotation, coupled by a torsional spring 
on the axis. This model is illustrated in figure 2. The engine-rotor 
inertia is represented as I e and the inertia of the articulated blades 

about the rotor shaft is represented as I r . The spring has a stiffness 

composed of two terms: the torsional stiffness of the drive shaft and 

a stiffness corresponding to the centrifugal force that tends to restore 
the rotor blades to radial positions. 

Consider first the balance of torques on the engine in figure 2: 


The aerodynamic torque Q e delivered by the engine to the rotating 
system is 


Simplified Model of Rotor 


Qe + Qs + Qp = 'iNe/cl't 


( 1 ) 


Q e = kyW + k e N e 

The torque Q s of the equivalent spring is 


( 2 ) 



( 3 ) 
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where 


and 




e = P (N e - N r ) dt 

The torque Q D developed by the damper is 

% = -k D (N e - N r ) = -k D de/dt 

Substituting the values of Q e , Q s , and Qp into equation (l) 


gives 


I e dN e /dt = k w W + k e N e + k s 9 - k D de/dt 


(4) 

(5) 

( 6 ) 

(?) 


The torque equation for the rotor is similarly derived. The aero- 
dynamic torque on the rotor is 

Or = k p p + k^j. (8) 

The rotor-torque equation becomes 

I r dN r /dt = kpS + k^ - k g e + k D de/dt (9) 

Equations (7) and (9) can be reduced to a convenient dimensionless 
form by dividing by The resultant equation can be expressed in 

terms of time constants and T r . Equation (7) becomes 


N, 


e p 


I r k W 1 


co r o) r I e (-k r ) 


W T r 1 ^e h Q _ Ir 




uu 


T e T 


jOij. ao r 


I e ( _ kr) 


e -2- 

CD- 


and equation (9) becomes 


N 


P 


N t 


S D 


1- = — P -L- -z- - -i— — + e + m a — 

oo r ^r “^r ^r ^i^r ^r l-kr; ^i^r u>r 


e -2- 


( 10 ) 


(ii) 


The reduction of equations (7) and (9) to dimensionless form has 
reduced the number of coefficients to be investigated. In equations (7) 
and (9) ; the coefficients are 


le> k W-> k e> k s-> k D> k p> k r> ^r 
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whereas in the dimensionless equation the coefficients are 

Ir i T*. *D % 

V W V V V k r 

Thus^ the number of coefficients have been reduced to six. The effects 
of only the first three are explored in this paper. 

Equations (5), (7), and (9) define the motion of the engine-rotor 
system as functions of fuel-flow disturbances W and rotor-pitch dis- 
turbance p. The coefficient k w is the slope of the torque - fuel- 

flow characteristic of the engine at constant speed. Similarly, k e is 

the slope of the engine-torque - speed characteristic at constant fuel 
flow, kp is the slope of the rotor-torque - pitch characteristic at 

constant speed, and k r is the slope of the rotor- torque - speed char- 
acteristic at constant pitch. The system of equations (5), (7), and (9) 
therefore constitutes a linearization of the motion of the engine and 
rotor. Such a linearized representation is applicable to small excur- 
sions from a steady- state condition and is useful for study of control 
stability. 

The system of equations (5), (7), and (9) or equations (5), (10), 
and (ll) can be simulated directly by electronic analog and therefore 
represents the extent to which the equations need be developed. Solution 
of these equations to give transfer functions among the inputs and out- 
puts will, however, yield useful information, because certain properties 
of the motion can be revealed by inspection of the equations . 


p 


Transfer Functions 

The equations (5), (7), and (9) can be solved to yield an equation 
showing the response of N r to W. If the derivative symbol d/dt is 

represented as p, this equation appears as follows: 


N 


r 


P k D - k s 

p^C + p2p + pE + G 


k W W 


( 12 ) 


when C, D, E, and G are constants having the following values: 

C = I e I r 

D = -I e k r - Ijvkg + I e kj) + IjJsp 

I e k s " k e k D “ k r k D 


E - -I r k s + k e k s 


( 13 ) 
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The remaining five transfer functions have the same 
different numerators. These transfer functions are: 

denominator, 

p 2 I e + p(k D - k e ) - k g 
Wr = p 3 C + p 2 D + pE + G kpP 

(14) 

p 2 I r + p(k D - k ) - k 
e p 3 C + p 2 D + pE + G w 

(15) 

pk D - k s 

N e “ p 3 C + p 2p + P E + G kpP 

(16) 

pl r - k 

0 = -= | kwW 

p°C + p^D + pE + G 

(17) 

-pl e + k e 

' p 3 C + p 2 D + pE + G kpP 

(18) 


Inspection of the denominator will, among other things, reveal the 
undamped natural frequency cn n and the pendulum frequency o) r of the 

rotor. If the energy dissipating coefficients k e , k r , and k^ equal 
zero, the denominator becomes: 

P 3 d e l r ) + P 2 (0) - p(l r k s + I e k s ) + 0 (19) 

Equating this term to zero gives the characteristic equation: 

P 2(l e I r ) - (I e + I r )k B = 0 (20) 

The undamped natural frequency is 

= Vfe +i )fe <2i) 

If the engine inertia I e is increased to infinity, the undamped 
natural frequency becomes the pendulum frequency of the rotor in 
torsion: 


CJD-r = 


(22) 
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The steady- state gains between inputs and outputs are readily 
obtained from the transfer function by evaluating the transfer function 
at zero frequency; this is done by factoring the numerators and denomi- 
nators so that their last terms are equal to unity. An important gain 
used in this report is the gain of engine speed to disturbance in fuel 
flow: 



(23) 


This gain is needed in analysis of control systems involving speed and 
fuel flow. 


SOLUTION OF EQUATIONS OF MOTION 
Analog Representation 

Solutions of the equations of motion were made with the help of an 
electronic analog. In this section, the analog representation used is 
presented, the manner of extending the representation to free-wheeling 
turbines is discussed, and the physical constants in the equation of 
motion are evaluated. 

The analog representation chosen for investigation of control is 
a direct representation of the elementary equations (5), (7), and (9) 
rather than a construction of the transfer functions (12), (14), (15), 
(16), (17), and (18). The method chosen is the simpler of these two, 
because it requires only three integrators (equaling the number of 
energy- storage elements), and it requires that each physical constant 
be represented only once. 

A direct representation of the equations of motion is shown in 
figure 3(a). In this figure, the transfer functions of the boxes are 
indicated. For example, the box containing kp is an amplifier with 

a gain numerically equal to kp, and the block with l/l e P is an inte- 
grator with an integrator time constant numerically equal to I e . It 

may be seen that the top row of boxes represents the equation of motion 
involving the engine speed N e j and the bottom row, the rotor speed N r . 

The subtracter and the box marked k s /p represent equation (5), and the 
remaining box, marked kp, represents the damper. 

The representation of a physical system in an electronic analog is 
usually made with all physical constants reduced to dimensionless form. 
Selection of scale factors is simplified, and the number of variables 
to be studied is reduced. In this study, the time scale was canceled 
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by multiplying by the rotor pendulum frequency co r . The resulting ana- 
log representation, shown in figure 3(b), was used in studies subse- 
quently presented in this report. In many applications, other reduc- 
tions to dimensionless form may be preferable. 

When the dimensionless representation is used in a control study, 
the loop gain of the control in terms of the dimensionless parameters 
is needed. The loop gain of the speed - fuel-flow control is consti- 
tuted as follows: 


In equation (24), the term Kq, which is the gain of the control, 

has units of the ratio of fuel-flow change to speed change. Equation 
(24) is produced by reducing the engine gain (eq. (23)) to dimensionless 
form and multiplying by Kq. 

In analyzing a control system, the second term on the right side of 
equation (24) is usually represented by an amplifier. The first term is 
already represented in figure 3(b). 


The analysis and the analog representation are easily extended to 
the case of the free-wheeling turbine drive. In this system, a gas- 
turbine engine provides hot gas at high pressure to an independent drive 
turbine that is geared to the rotor shaft. The operating conditions of 
the gas-turbine engine are approximately independent of the load applied 
to the drive turbine. 

The analog representation of the rotor and its free-wheeling drive 
turbine is made by a simple replacement of the fuel-flow input by the 
two new inputs, gas pressure and temperature to the drive turbine. The 
analog representation of the response of gas pressure and temperature 
to fuel flow supplied to the gas turbine is taken from figure 4(b) of 
reference 4. 


1 



Extension to Free-Wheeling Turbine Drives 


Evaluation of Physical Constants 


Torsional coefficient k s . - 


Coefficients in equations (5), (7), 
evaluation are k s and k r . In 


and (9) that pose difficulties in 
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particular, the torsional elasticity k s is important, because the ade- 
quacy of the assumption of pure rotation depends upon the accuracy with 
which k s can be chosen. 

The total torsional deflection amplitude of the center of gravity 
of a blade, relative to the engine rotor, is 6. This deflection is 
composed as follows: 


DO 

CD 


® a (25) 

where £ is the amplitude of angular motion of the blades and a is 
the amplitude of twist of the drive shaft. Now 

i = Qs/ k a (26) 

» = QsAb (27) 


therefore, 



In figure 2, k s was defined by the relation 


(28) 


k s 



(29) 


Combining equations (28) and (29) and eliminating Q s /e yield 



(30) 


The torsional elasticity k^ of the shaft is easily determined) 

evaluation of the rotor deflection constant k a remains to be 
accomplished. 


The geometry of the force vectors acting on a rotor blade is shown 
in figure 4. The rotor blade is hinged to the hub at a distance a 
from the center of the drive shaft. In steady rotation, the centrifugal 
force F r acting through the center of gravity of the blade attempts to 

hold the blade so that its center of gravity is on a line with the hinge 
and the center of the rotor shaft. An aerodynamic drag force Fjj nor- 
mal to the radius causes the blade to swing back to an angle £ where 
the tensile force Ft along the axis of the blade has a component par- 
allel, equal, and opposite to the drag force. 
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The definition of k a is 

= d ^sM (31) 

Now 

«s = F t f * 5 F r f 
F r = M(d + e) 

and 


f = a sin £ 


For small values of 


f = a 5 


Hence 


=t2 


Q s = MN r (d + e) £ 


whence 


dQ £ 

W 


= K 


■ ( 4 ^)(!) 


a 2 MN^ 


A close approximation having more convenient dimensions is 

k a - (H^) “ 2 ™r ( 32 > 

Damping coefficient k r . - The rotor damping coefficient k r is 

easily found if the shaft horsepower in steady state is known. The 
definition of k r is 


k r - 



(33) 


where the subscript indicates that the rotor pitch P is held constant 
while the slope of the torque-speed curve is measured. If consideration 
is given to the fact that at constant rotor pitch the torque is propor- 
tional to the square of the rotor speed. 


k r 


-2X550 

N r 


(34) 
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In the dimensionless representation of figure 3(b), the rotor time 
constant t r is used. The rotor time constant is 

I r _ I r N^ 

Tr " -}c r " 2X550 hp ( 35 

COMPARISON OF CALCULATED AND EXPERIMENTALLY 
OBSERVED TORSIONAL MOTION 

A critical comparison of torsional motion calculated by use of 
equations (5), (7), (9), and (32) with experimental observation is 
necessary, because these equations are based upon approximations. 
Experimental data of the type needed for this comparison are contained 
in reference 2. 

In reference 2, measurements of the torsional oscillatory response 
of a helicopter rotor to a sinusoidal disturbance in rotor pitch are 
reported. The tests were performed on a rotor from a 2500-pound heli- 
copter . The rotor was mounted on a helicopter test tower and driven by 
a 1350-horsepower reciprocating engine. Pertinent details of the rotor 
dimensions and test equipment are given in table I. The data used for 
comparison in the present report are the data for no oil damping shown 
in figure 2 of reference 2. 


Calculation of Rotor Motion 

The analog representation of figure 3(b) was used in calculating 
the rotor motion. The following coefficients were determined using 
the data from table I: 


1 Coefficient 

Value 

k a 

-22,800 (ft- lb) (sec) /radian 

k s 

-15,000 ft-lb/ radian 

k r 

-263 (ft-lb) (sec) /radian 

(JD r 

5.22 radians/sec 

T r 

2.10 sec 

1/ (a r a>r) 

0.0913 

I r/ I e 

1.34 

k D 

0 (ft-lb) (sec) /radian 
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The engine damping term k e (or t 0 ) is not reported in reference 2j 

it therefore had to he estimated. This estimation was made by assuming 
the indicated horsepower at constant throttle setting to be independent 
of engine speed and the friction torque of the engine to be proportional 
to the square of the engine speed f In the test, the engine, which was 
operated at 55 percent of the speed for maximum power, delivered 9.6 
percent of its maximum power. The friction power was computed for these 
conditions on the assumption that the friction was in the same propor- 
tion as that in the Plymouth engine shown in figure 311, page 49, of 
reference 5. The computation yielded the following results: 

k r = - 756 (ft-lb) sec/radian 
T e = 0.542 sec 


Hence, 


I r Terror) 


0.26 


Comparison with Experimental Observation 

The calculated and the observed torsional motion of the rotor are 
compared in figure 5 where drag-angle amplification is shown as a func- 
tion of the frequency of the sinusoidal disturbance in rotor pitch. 

The drag-angle amplification is defined as the ratio of the ampli- 
tude of drag-angle oscillation £ at the selected disturbance fre- 
quency a) to the change in the steady- state value of the drag angle 

at constant rotor speed when a change in rotor pitch of the same 

•magnitude as the amplitude of the sinusoidal disturbance in pitch is 
made. The rotor-pitch frequency cd is made dimensionless by dividing 
by Np, which is constant at 23 radians per second. The circles in- 
dicate data points reported in reference 2, and the curve shows the 
response computed using the analog representation of figure 3(b). 

The data are in excellent agreement with the computed response. The 
resonant frequency is in close agreement - 0.35 times the rotor speed. 

The damping, as indicated by the drag-angle amplification factor at res- 
onance, shows much better agreement than can normally be expected in 
view of the uncertainties in estimating the damping terms. 

A disagreement between calculated and observed response occurs at 
high frequency, above a dimensionless frequency of 0.42. The increasing 
divergence between calculated and observed response at higher frequen- 
cies is interpreted as an indication that a system of differential 


2977 


2977 


NACA TN 3027 


15 


equations of higher than the third order represented in equations (5), 
(7), and (9) is needed, because the rotor blades can oscillate about 
their own centers of gravity in addition to having their centers of 
gravity oscillate around the rotor shaft. 

In general, the agreement between calculated and observed response 
indicates that the differential equations (5), (7), and (9) and the ana 
log representations in figure 3 will give satisfactory results in the 
analysis of engine control. With very large helicopters (50,000 lb 
gross wt) having heavy rotor blades, the resonant frequency will be so 
low that the error above the resonant frequency may become serious . 
Further investigation of the representation of the frequency response 
above resonance and additional careful measurement of response in this 
frequency range are needed. 


Other Responses of the Rotor of Reference 2 

As a matter of interest, other responses of the rotor described in 
reference 2 were computed. The frequency responses of engine speed and 
rotor speed (blade centers of gravity) to sinusoidal changes in rotor 
pitch and the transient responses to a step change in rotor pitch were 
calculated on the electronic analog. These responses are presented in 
figure 6 in the form of photographs of an oscilloscope screen. 

The frequency response to rotor pitch disturbance is shown in fig- 
ure 6(a). The blade-angle input in figure 3(b) was excited by a sweep- 
frequency oscillator that produced a sine wave, the frequency of which 
increased linearly with time. The x-axis of the oscilloscope was so 
synchronized with the sweep-frequency generator that, in the photograph, 
the x-axis is proportional to frequency. The y-axis was attached to the 
proper output, thereby creating an oscilloscope image comparable with 
the frequency response in figure 5, with the exception that the scales 
are linear rather than logarithmic. The drag-angle response shown 
should be the same as previously presented in figure 5. 

The frequency response of engine speed is interesting in that a 
pronounced minimum occurs between zero frequency and resonance. At 
this "antiresonance" frequency, disturbances in rotor pitch have prac- 
tically no effect on engine speed. 

The corresponding transient response to a step change in rotor 
pitch is shown in figure 6(b). The highly underdamped oscillatory 
motion superimposed upon the exponential drift toward equilibrium is 
obvious. The oscillations in torque have an amplitude almost equal to 
the total changes in torque produced. 
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CONTROL STABILITY 

The stability of engine speed and torque controls was investigated 
using the analog of figure 3(b) as a representation of the engine. 
Before making this study, however, it was necessary to estimate the 
proper dynamic coefficients for a large rotor to match a gas-turbine 
engine for which the dynamic coefficients were known. In this section 
of the report, the effects of helicopter weight on rotor dynamics are 
first investigated, and a rotor to match a 2700-horsepower gas turbine 
is chosen. Then, the stability of engine-speed and torque controls are 
explored. 


Effect of Helicopter Weight on Rotor Dynamics 

The gas-turbine engine for which the dynamic coefficients are known 
produces 2700 horsepower. The rotor for which the dynamic properties 
are known (the rotor studied in ref. 2) should be used with a 210- 
horsepower engine. This rotor was tested with a 1350-horsepower recip- 
rocating engine. If the power-weight ratio of the helicopter is assumed 
to be 0.085 horsepower per pound regardless of weight, the 2700- 
horsepower turbine would be suitable for a 32,000-pound helicopter. 

The dimensions and properties of a rotor suitable for a 32,000- 
pound helicopter may be computed by use of the data presented in table 
II. In this table, two different laws of rotor weight are used: One 

law makes the blade weight proportional to the helicopter plan form 
area, whereas the other law makes the blade weight proportional to the 
three-halves power of the helicopter weight. The first law is probably 
closer to what may be accomplished, but the results are somewhat opti- 
mistic in that the weight of a large rotor would be underestimated. 

The effect of helicopter gross weight on the rotor dynamic coef- 
ficients is shown in figure 7. This figure was computed using the 
dimensions of the rotor for the 2500-pound helicopter and assuming the 
law of direct proportionality in table II. The dimensions of a rotor 
for a 32,000-pound helicopter were estimated, and the characteristics 
of a gas-turbine engine having a time constant of 1 second and an iner- 
tia of 1,15 pound-foot second^ (applicable to 2700 hp engine) were 
assumed. The engine time constant and the ratio of the inertia of the 
engine to that of the rotor were then assumed to be independent of heli- 
copter gross weight. The resonant frequency and drag angle and the fre- 
quency response of engine speed to fuel flow were determined with the aid 
of the electronic analog of figure 3(b). The results presented in fig- 
ure 7 show that the drag-angle amplification, engine-speed amplifica- 
tion, and resonant frequency decrease with increasing gross weight of 
the helicopter. The significant parameters for a gross weight of 
32,000 pounds are listed in the following table together with those 
appli< >le to the test equipment used in securing the data of figure 5: 
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Parameter 

Helicopter gross weight, lb 


32,000 

2500 

Radius of rotor, ft 

68 

19 

Rotor speed, radians/sec 

6.43 

23.0 

a) r , radians/sec 

1.46 

5.22 

Engine inertia at engine speed. 

1.15 


lb-ft sec 2 



Engine speed, radians/sec 

1497 


Engine inertia at rotor speed I e , 

62,000 


lb-ft sec 2 



Ir/Ie 

1.45 

1.34 

l/ (t ru>r ) 

0.324 

0.0914 

I e T r 



I r 

0.462 

0.26 


Comparison of the properties of the large helicopter and the small 
helicopter cannot be made because the engine for the smaller helicopter 
was a reciprocating engine and that for the large helicopter, a turbine- 
propeller engine. 

Control system . - The control systems studied are shown in figure 
8 - Figure 8(a) shows a proportional -plus-integral control of engine 
speed in which speed is controlled by manipulation of fuel flow. An 
integrator time constant of 1.45 seconds and a fuel-system lag 

t c of 0.2 second were used throughout the studies. The rotor system 

was simulated by the analog of figure 3(b). Except where otherwise 
noted, the properties of the 32,000 pound helicopter were used. 

Effect of torsional motion on control response . - A direct impres- 
sion of the effect of the torsional motion on control response and sta- 
bility can be obtained by comparing the responses of an articulated 
rotor and the response of a rigid rotor system having the same inertia 
and power input. This comparison is made in figure 9. Figure 9(a) 
shows the transient response of the articulated rotor and figure 9(b), 
the response of the rigid rotor. Computations for the rigid rotor are 
given in the appendix. The control used in each case had a loop gain 
of 2.5. 

The oscillatory contribution of the articulated rotor is evident 
in a comparison of figures 9(a) and 9(b) . A much larger overshoot in 
torque occurs with the articulated rotor than with the rigid rotor, 
and the time to establish a new equilibrium is much longer with the 
former system. 
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As a matter of interest, the transient responses of the engine 
rotor systems with articulated and rigid rotors without engine controls 
are shown in figures 9(c) and 9(d). The articulated rotor exhibits an 
underdamped high-frequency oscillation superimposed upon an exponential 
response, whereas the rigid rotor exhibits only a smooth exponential 
response . 

It is concluded that stable control of engine speed by fuel flow 
can be established, but that the control system will be more oscillatory 
than it would be if used on a perfectly rigid rotor system. 

Selection of engine -speed control . - The control with the high 
loop gain, 2.5, needed for accurate speed control shows a 22 percent 
overshoot in torque. The overshoot is determined by scaling the torque 
in figure 9(a). This torque overshoot may be practically eliminated by 
reducing the gain, as shown in figure 10. In this case the loop gain 
was reduced to 0.9. The low loop gain is shown to have the disadvantage 
of permitting a large excursion in engine speed. The speed excursion at 
the low gain is two times the speed excursion at high gain. 

Effect of gross weight on speed control . - The effects of helicop- 
ter gross weight on dynamic response as well as the dynamic coefficients 
are shown in figure 7 . It may be seen that not only did the resonant 
frequency decrease with increasing gross weight (as would be expected), 
but also that the damping of the resonant oscillation increased, as 
revealed by the decreasing trend in the ratio of resonant amplification 
to steady-state amplification. This effect is a result of the assump- 
tion that the rotor-blade weight is proportional to blade plan form 
area . 


The corresponding variation of control behavior with change in 
gross weight is shown in figure 11. The loop gain above which continu- 
ous oscillation occurs is shown to increase to a value of 3.7 at a 
gross weight of 20,000 pounds, and further increase in gross weight 
has no effect on this stability limit. This decrease in loop gain at 
the low gross weight is attributable to the comparatively low damping 
of the rotor. With the high gross weight, the loop gain is only 3.7j 
whereas, for a rigid rotor, the loop gain is 17. Thus, the torsional 
motion of the rotor contributes greatly to the instability of the engine 
control. 

Effect of rotor blade weight on control . - As pointed out in the 
section on the Effect of Helicopter Weight on Rotor Dynamics, the 
method used to estimate rotor weight might give low estimated rotor 
weights for large helicopters. In a method that gives weights that 
are probably too high for large helicopters, the weight of the rotor 
is assumed to be proportional to the volume of the rotor blade and, 
hence, proportional to the three-halves power of the helicopter gross 
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weight. Formulas for the dynamic constants under this assumption are 
given in table II. Extrapolation of the inertia of the small rotor of 
reference 2 to a size adequate for a 32,000-pound helicopter results in 
an inertia I r /I e of 51.8 as con^ared with 1.45 when the law of direct 
proportion is assumed. 

The effect of the inertia ratio on the maximum stable loop gain of 
the engine -speed control is shown in figure 12; the dimensions of the . 
rotor were held constant, and only the weight and inertia were varied. 
As the rotor weight is allowed to increase, the maximum permissible 
loop gain decreases. Analysis for inertia ratios above 11.6 was impos- 
sible with the particular analog equipment used. 

Transient responses of a helicopter having a rotor four times the 
mass (inertia ratio of 5.8) of that used in figure 9 are shown in fig- 
ure 13. Figure 13(a) shows that the torsional oscillation is much less 
damped than that of the lighter rotor ; the speed responses are, of 
■Q course, much slower than with the lighter rotor. 

■o 

| Tor< 3 ue control . - The torque control shown in figure 8(b) was in- 

vestigated. This control had the same components and coefficients as 
the speed control except that it was found necessary to add derivative 
action to the calculator. The derivative time constant T d was ad- 
justed to meet the stability requirements of the control. 

The torque control proved to be very prone to oscillate. Diver- 
gent oscillatory instability occurred above a loop gain of 2.9 and at 
lower loop gain the system showed oscillation. 

A very smooth response was obtained at the high loop gain of 4.1 
when a derivative time constant of 0.6 second was used. The system 

was very sensitive to both gain and derivative time constant and would 
break into continuous oscillation with either a higher* or a lower de- 
rivative time constant. 

Although a torque control could be effective, it would have a 
strong tendency to oscillate. 


CONCLUSIONS 

Analysis was made of the effect of the torsional oscillations in 
hinged rotor blades in a helicopter on the stability of engine controls 
on a gas-turbine engine geared to the rotor. The analysis led to the 
following conclusions: 
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1. Torsional oscillation of the articulated blades contribute con- 
siderably to the instability of engine controls. The unstabilizing 
action is more evident in torque control by fuel flow than in engine- 
speed control by fuel flow. 

2. Increasing the mass of rotor blades reduces the stability of 
control systems, thereby necessitating lower loop gains and resulting 
in greater excursions in rotor speed. 

3. The rotor-engine system may be represented adequately by two 
rotating masses connected by a torsional spring. The resonant fre- 
quency and damping were accurately represented, but the amplitude re- 
sponse above resonance was underestimated. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 12, 1953 
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APPENDIX - RESPONSE OF RIGID ROTOR 


The response of a rigid rotor can be derived by setting the deflec- 
tion 9 in equations (5), (7), and (9) equal to zero and solving for 
N e . The resultant equations are 

N e ( I e + J r)P = k pP + k w w + ( k e + k r ) N e ( A1 ) 

-Q s = k w W + k e N e - I e N e p (A2) 

These equations, when reduced to the same dimensionless parameters 
used for the articulated rotor, become: 




k w 

W 1 

k p p 

N e 1 

v r . ii 

CO r 0) r l_I r J 
tr\ '\ 

1 T r CD r 

TvT 

1 

oo r T I ^> r 

(-k r ) cn r 

0) r T 1 0> r 

LVe J 

(A3) 

/V 

Vrj 

1 kw 

i 

W leTr 

1 N e 

Ie Ne p 

(A4) 

<o| 

(-k r ) 

XjCDp CDp IpX 0 


Ir 


Analog representations of these equations are shown in figure 14. 
Figure 14(a) is a direct representation in real time and is a construc- 
tion of equations (Al) and (A2). Figure 14(b) is a nondimensional rep- 
resentation and is a construction of equations (A3) and (A4). 

For the 32,000-pound helicopter, the following coefficients were 
found: 


I e 

Y~ + 1 = 1.69 
1 r 



= 0.795 



= 0.472 
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TABLE I . - DIMENSIONS AND PROPERTIES OF HELICOPTER ROTOR AND 


ENGINE USED IN TEST 
jpata from ref . 2 



Helicopter rotor: 

Blade radius, ft 

Blade weight, 3 blades, including drag-hinge assembly, slugs . 7.7 

Blade moments of inertia about their centers of gravity 

(3 blades), lb-ft sec 2 . 2 24 

Biade center-of -gravity location from vertical pin axis, ft . . 5.76 
Vertical -pin-axis location from center of rotation, ft . . 0 757 

Drive-shaft torsional stiffness, ft-lb/radian .... 43’eoo 

Effective moment of inertia of engine and gearing refer- 
enced to rotor shaft speed, lb-ft sec 2 4 10 

Rotor shaft speed, radians/sec p. 

Shaft horsepower 9 


Engine: 

Rated horsepower . 


1350 

2400 

1320 

2500 


Rated engine speed, rpm 


Engine speed at test conditions 
Helicopter gross weight, lb 


, rpm . 
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TABLE II. - EFFECT OF HELICOPTER GROSS WEIGHT ON ROTOR 

DIMENSIONS AND PROPERTIES 


Coefficient 

Proportional 

(a) 

Three-halves 

(b) 

Horsepower 

C L 

L 

I-r 

L 2 

L 5 / 2 

Jl 

r < 

L 1 / 2 

Ll/2 

t 

k Q 

L 

L 3 / 2 

xv a 

L 

L 3 / 2 

s 

CD r 

i/l 1 / 2 

1/LV2 

-L 

ky. 

L 2 

L 2 

r 

T 

Constant 

L 1 / 2 

r 

*e 

L 2 

L 2 

Tp 

Constant 

Constant 

I r/ I e 

Constant 

L 3 / 2 

l/ (t^) 

L 1 / 2 

Constant 

1 e T r 

X r ^e^r^r) 

L 1 / 2 

l/L 

XeAr 

Constant 

l/L 1 / 2 


a Rotor weight assumed proportional to plan 
form area. 

^Rotor weight assumed proportional to three- 
halves power of blade plan form area. 
c Gross weight of helicopter. 
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DO 
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Engine Rotor 

S / N 



where 9 = / (N e - N r ) dt 

and Q s = k s /(N e - N r ) dt = k s © 



Figure 2. - Representation of engine and rotor as simple system of bodies in pure 
rotation. 
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(a) Direct representation of system in real time. 

Figure 3. - Analog representation of articulated helicopter rotor and engine. 
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(t>) Nondimensional representation. 



Figure 3. - Concluded. Analog representation of articulated helicopter rotor and engine 
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Figure 5. - Frequency response of blade drag angle to disturbance in rotor- 
pitch angle. 


2977 


2977 


Engine 

speed 


Eotor 

speed 


Drag angle 
or torque 


0 



Torque 


Engine 

speed 


Eotor 

speed 


.568 
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7.7 


Eotor pitch frequency, o) / N. 

(a) Frequency response to disturbance in 
rotor pitch. 


Time, sec. 

(b) Transient response to step disturbance 
in rotor Ditch. 


Figure 6. - Example of results obtained with analog representation of rotor-engine system presented in figure 2 


OJ 

H 


NACA TN 3027 




Resonant amplification e T r 1 

Frequency, radians/sec Steady-state amplification I r T e (t r ox r ) 'tjXUp 
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Gross weight of helicopter, lh 


Figure 7. - Effect of helicopter weight on dynamic characteristics of rotor. 
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(a) Engine-speed control by fuel flow. 



Figure 8. - 


(b) Engine-torque control by fuel flow. 
Control systems investigated ( T . , 1.45 sec; T , 



0.2 sec). 
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(a) Articulated, rotor with engine control. 


0 Time, sec 27 

(1) Rigid rotor with engine control. 


Figure 9. 
pitch . 


- Comparison of transient responses of an articulated and a rigid helicopter rotor to step change in rotor 
Helicopter weight, 32,000 pounds; control- loop gain, 2.5. 
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Time, see 
(c) Articulated rotor, no control. 
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Time, sec 

(d) Rigid rotor, no control. 


Figure 9. - Concluded. Comparison of transient responses of an articulated and a rigid helicopter rotor to step change 
in rotor pitch. Helicopter weight-, 32,000 pounds; control- loop gain, 2.5. 
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Figure 10. - Transient response of articulated rotor system with control- loop gain 
adjusted to highest giving no significant overshoot in torque. Response to step 
change in rotor pitch; control- loop gain, 0.9. 


2977 


Control loop gain above which 

continuous oscillation occurs 


NACA TN 3027 


37 


CD 

00 



Figure 11. - Effect of helicopter gross weight on stability 
of engine -speed controlled by fuel flow. 


Maximum stable loop gain 



Figure 12. - Effect of rotor weight (inertia) on stability limit of engine-speed 
control in helicopter weighing 32,000 pounds. Control, integral plus propor- 
tional; integrator time constant, 1.45; control lag, 0.2 second. 
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(a) No control. 
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(b) Control with maximum stable loop 
gain of 0.7. 


Figure 13. - Transient responses of 32,00-pound helicopter having rotor with four times -weight represented in figure S. 
Inertia ratio, 11.6. e 
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(b) Nondimensional representation. 

Figure 14. - Analog representation of rigid helicopter rotor and engine. 
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